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Abstract 

We present the first measurement of the D*^ width using 9/fb of e^e~ data 
collected near the T(4S) resonance by the CLEO II. V detector. Our method 
uses advanced tracking techniques and a reconstruction method that takes 
advantage of the small vertical size of the CESR beam spot to measure the 
energy release distribution from the D*^ D^tt'^ decay. We find T{D*^) = 
96 lb 4 (Statistical) it 22 (Systematic) keV. We also measure the energy release 
in the decay and compute Am = m£,*+ —niDO = 145.412ib0.002(Statistical) it 
0.012 (Systematic) MeV/c^. 
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I. INTRODUCTION 



A measurement of r{D*~^) opens an important window on the non-perturbative strong 
physics involving heavy quarks. The basic framework of the theory is well understood, 
however, there is still much speculation - predictions for the width range from 15keV to 
150 keV The level splitting in the B sector is not large enough to allow real strong 
transitions. Therefore, a measurement of the width of the D*~^ gives unique information 
about the strong coupling constant in heavy-light meson systems. 

The total width of the D*~^ is the sum of the partial widths of the strong decays D*^ —>■ 
D°7r+ and D*^ D^n^ and the electromagnetic decay D*^ D^^. We can write the width 
in terms of strong couplings, gD*^D"-K+ and go'^D+wO, and an electromagnetic coupling, 

gD*^D+'y- 

r(D*+) = T{D\+) + r(Z}+7r°) + r(D+7) (1) 

2 2 2 

where the momenta are those for the indicated particle in the D*~^ rest frame, and a is the 
fine structure constant. This can be rewritten using the isospin relationship 



gD*^D-n 



— ^go-'-^D+TTO = gD*^DOTT+, (3) 



and relating gD*~*D-K to a universal strong coupling between heavy vector and pseudoscaler 
mesons to the pion, g, with 

gD*^D-K = — g, (4) 
where fir is the pion decay constant. All this yields 



o 2 2 2 

r(n*+\ - ^3 3 ,_g 3 I ^9d^^d+^ 3 , . 



The width of the D*^ only depends on (7 p| since the contribution of the electromagnetic 
decay with branching fraction (1.68 ± 0.45)% [Q can be neglected. The measurement of g is 
needed in the extraction of Vub in semileptonic b ^ u decays . 

Prior to this measurement, the D*~^ width was limited to be less than 131 keV at the 
90% confidence level by the ACCMOR collaboration ||^. The limit was based on 110 signal 
events reconstructed in two decay channels with a background of 15%. This contribution 
describes a measurement of the D*^ width with the CLEO II. V detector where the signal, in 
excess of 11,000 events, is reconstructed through a single, well-measured sequence, D*^ 
7r^Q^D°, —>■ K~7T~^. Consideration of charge conjugated modes are implied throughout 
this paper. The level of background under the signal is less than 3% in our loosest selection. 

The challenge of measuring the width of the D*~^ is understanding the tracking system 
response function since the experimental resolution exceeds the width we are trying to mea- 
sure. Candidates with tracks that have mismeasured hits, errors in pattern recognition, and 
large angle Coloumb scattering are particularly dangerous because the signal shape they 



4 



project is broad and the errors for these events can be underestimated, resulting in events 
that can easily influence the parameters of a Breit-Wigner fitting shape. We generically term 
such effects "tracking mishaps." A difficulty is that there is no physical calibration for this 
measurement. The ideal calibration mode would have a large cross-section, a width of zero, 
decay with a rather small energy release to three charged particles one of which has a much 
softer momentum distribution than the other two which decay through a nearly zero width 
resonance with a measurable flight distance. Such a mode would allow us to disentangle 
detector effects from the underlying width but no such mode exists. 

Therefore, to measure the width of the D*^ we depend on exhaustive comparisons be- 
tween a GEANT p based detector simulation and our data. We addressed the problem by 
selecting samples of candidate D*^ decays using three strategies. 

First we produced the largest sample from data and simulation by imposing only basic 
tracking consistency requirements. We call this the nominal sample. 

Second we refine the nominal sample selecting candidates with the best measured tracks 
by making very tight cuts on tracking parameters. There is special emphasis on choosing 
those tracks that are well measured in our silicon vertex detector. This reduces our nominal 
sample by a factor of thirty and, according to our simulation, has negligible contribution 
from tracking mishaps. We call this the tracking selected sample. 

A third alternative is to select our data on specific kinematic properties of the D*^ decay 
that minimize the dependence of the width of the D*^ on detector mismeasurements. The 
nominal sample size is reduced by a factor of three and a half and, again according to our 
simulation, the effect of tracking problems is reduced to negligible levels. We call this the 
kinematic selected sample. 

In all three samples the width is extracted with an unbinned maximum likelihood fit 
to the energy release distribution and compared with the simulation's generated value to 
determine a bias which is then applied to the data. These three different approaches yield 
consistent values for the width of the D*^ giving us confidence that our simulation accurately 
models our data. 

II. CLEO DETECTOR AND DATA SAMPLES 

The CLEO detector has been described in detail elsewhere. All of the data used in 
this analysis are taken with the detector in its II. V configuration [0. This work mainly 
depends on the tracking system of the detector which consists of a three-layer, double sided 
silicon strip detector, an intermediate ten-layer drift chamber, and a large 51-layer helium- 
propane drift chamber. All three are in an axial magnetic field of 1.5 Tesla provided by a 
superconducting solenoid that contains the tracking region. The charged tracks are fit using 
a Kalman filter technique that takes into account energy loss as the tracks pass through the 
material of the beam pipe and detector 0. 

The data were taken in symmetric e+e^ collisions at a center of mass energy around 
10 GeV with an integrated luminosity of 9.0/fb provided by the Cornell Electron-positron 
Storage Ring (CESR). The nominal sample follows the selection of D*^ — > vr^^^D^ — > 
K^Ti^-n^Y^^ candidates used in our — mixing analysis |]^. 
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Our reconstruction method takes advantage of the small CESR beam spot and the kine- 
matics and topology of the D*~^ tt^^^D'^ — > 7i^^-^K~7i~^ decay chain. The K~ and vr"*" are 
required to form a common vertex. The resultant candidate momentum vector is then 
projected back to the CESR luminous region to determine the production point. The 
CESR luminous region has a Gaussian width ~ 10 fim vertically and ~ 300 fim horizontally. 



It is well determined by an independent method . This procedure determines an accurate 
D° production point for D^'s moving out of the horizontal plane; D°'s moving within 0.3 
radians of the horizontal plane are not considered. Then the vr^o„ track is refit constraining 
its trajectory to intersect the production point. This improves the resolution on the 
energy release, Q = M^K^n^n^i^^) — M^K^n^) — m^+, by more than 30% over simply 
forming the appropriate invariant masses of the tracks. The improvement to resolution is 
essential to our measurement of the width of the Our resolution is shown in Figure |I] 

and is typically 150 keV. The good agreement between Monte Carlo and data demonstrates 
that the kinematics and sources of uncertainties on the tracks, such as the number of hits 
used and the effects of multiple scattering in detector material, are well modeled. 

To further improve the quality of reconstruction in our nominal sample, we apply some 
kinematic cuts to remove a small amount of misreconstructed signal and background. Fig- 
ure 1^ shows the distribution of the momentum of the vr^ow ^ function of the D*~^ candidate 
momentum. We apply a cut at the kinematic boundary as shown in the figure. Figure ^ 
also shows the opening angle 9 between the vr^^^ and the candidate as a function of the 
D*'^ candidate momentum. We apply a cut of < 38°/PD*[GeV] which is just beyond the 
kinematic limit to account for resolution smearing. We also require ag < 200 keV which 
removes the long tail in the error distribution. 

The tracking selected sample makes much more stringent cuts on the quality of the tracks 
used to identify the candidates. All tracks are required to have hits in both the r0 and z 
views in all three layers of the silicon strip detector as opposed to the nominal two silicon 
hits per view. None of these hits are allowed to be within 2 mm of a silicon wafer edge. The 

daughter tracks are required to have at least 38 of the possible 51 main drift chamber 
hits and seven of the ten intermediate drift chamber hits. The P^r degree of freedom of 
the fit to these two tracks are limited to less than 2 in each of the two drift chambers and 
50 in the silicon strip detector. These selections are designed to remove tracks that have 
tracking mishaps or decay in flight. 

We compare the simulation and the data as a function of kinematic variables of the D*^ 
decay. This will provide another test of the simulation's modeling of the data, and be the 
basis of our study of systematic uncertainties in the analysis. The most important kinematic 
variables are the "derivatives" which are defined by 



= M{K7Tf + ml+ +2{EnoE^+ -Pd»P^+ cosO) , (6) 

slow slow slow 

Pdo=Pdo/Edo, (7) 

slow slow slow 

dQ 



dPno M 



if3D0-f3,+ cos 6), (9) 
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FIG. 1. Distribution of ctq, the uncertainty on Q as determined from propagating track fitting 
errors. The arrow indicates a selection discussed in the text. 
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FIG. 2. The ir^i^^ momentum (a) and the opening angle between and vr^'^^^ (b) both versus 
the D*^ momentum in the nominal data sample. 



-(3nocose), (10) 



slow 



= ^sin6'. (11) 

These derivatives test correlations among the basic kinematic variables, the and tt^q^ 
momenta and the opening angle, 6. We compare by dividing the Q distribution into ten slices 
in each of the kinematic variables and fitting the ten sub-distributions of Q to Gaussians. 
We display the width and mean of the ten fits as a function of each of the six kinematic 
variables in Figures ^ and ^. 

The quality of the width comparison (Figure]^) is excellent, with the simulation generated 
with an underlying r(D*+) in the range of 90-100 keV agreeing well with the data for all 
the kinematic variables. Even when generated with an underlying r{D*~^) = keV the 
simulation accurately follows the data's changes as the kinematic variables vary across their 
allowed range. 

The quality of the mean comparison (Figure ^) is not as good. The dependence of the 
mean of Q is not well modeled versus the n^ow momentum, dQ/dP^+ , and dQ/dPoo by 

slow 

our simulation. We discuss the consequences of this imperfect modeling of the data in the 
section on systematic uncertainties below. 

Figure ^ shows the three derivatives plotted against each other in the data. Note that 
if we select dQ/dP^o and dQ/dP^+ both to be close to zero we minimize the dependence 

slow 

of Q on the basic kinematic variables Pdo and P^+ , and thus minimize the contribution of 

slow 

the kinematic variables to the width of the Q distribution. With this selection we are more 
sensitive to the underlying width of the Q distribution rather than variations caused by any 
mismodeling of Q's dependence on the basic kinematics. The kinematic selection is defined 
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FIG. 3. Gaussian width of Q distribution versus kinematic parameters and derivatives. 
Data; o-Simulation with T^*^ = 0; A-Simulation with r/)*^ = 90 keV. 
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FIG. 4. Gaussian mean of Q distribution versus kinematic parameters and derivatives. 
•-Data; A-Simulation with Tj:)*^ = 90 keV. The horizontal Unes show the average value of Q 
for the two samples. 
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dQ / dRpO 

FIG. 5. Distributions of dQ/dPjjo versus dQ/dP^+ (a), dQ/dO versus dQ/dP^+ (b), and 

slow slow 

dQ/dO versus dQ/dP£,o (c) in the data. The dashed regions show the selection that defines the 
kinematic selected sample. 
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dQ 



< 0.005, 



< 0.05. 



(12) 
(13) 



Table |T| summarizes the statistics in our three samples. The tracking and kinematic 
samples are subsets of the nominal sample. The two subsets contain 94 common candidates. 



III. FIT DESCRIPTION 

We assume that the intrinsic width of the is negligible, r{D^) ^ r{D*^), implying 
that the width of Q is simply a convolution of the shape given by the D*~^ width and 
the tracking system response function. Thus we consider the pairs of Q and ctq for D*^ 
7r^Q^D° K'n^Ti^^^ where aq is given for each candidate by propagating the tracking errors 
in the kinematic fit of the charged tracks. We perform an unbinned maximum likelihood fit 
to the Q distribution. 

The underlying signal shape of the Q distribution is assumed to be given by a P-wave 
Breit-Wigner with central value of Q, Qq. We considered a relativistic and non-relativistic 
Breit-Wigner as a model of the underlying signal shape, and found negligible changes in the 
fit parameters between the two. The width of the signal Breit-Wigner depends on Q and is 
given by 

where Fq is equivalent to F (/}*+), P and M are the candidate vr^ow momentum in the 

D*^ rest frame and Kt^t^sIow mass, and Pq and Mq are the values computed using Qq. The 
effect of the mass term is negligible at our energy. The partial width and the total width 
differ negligibly in their dependence on Q for Q > 1 MeV. 

For each candidate the signal shape is convolved with a resolution Gaussian with width 
ctq, determined by the tracking errors, as a model of our finite resolution shown in Figure |I[ 

The fit also includes a background contribution with a fixed shape. The shape for the 
background is taken from fits to the background prediction of our simulation with a third 
order polynomial. The level of the background is allowed to float in our standard fit. The 
predicted background shape and fits are displayed in Figure ^. 

Figure 1^ shows the Q distribution for our nominal data sample. Note that besides the well 
measured signal and the small, slowly varying background, there is also a small component 
centered on the signal with a large width. Therefore we allow a small fraction of the signal, 
fmis, to be parametrized by a single Gaussian resolution function of width amis- This shape 
is included in the fit to model the tracking mishaps which our simulation predicts to be at the 
5% level in the nominal sample and negligible in both the tracking and kinematic selected 
samples. Typically we constrain the level of this contribution while allowing amis to float. 
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FIG. 6. Our simulation's prediction of the background for the three samples discussed in the 
text. Also shown are the fits to third order polynomials that are used in the fits to the data. 
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FIG. 7. The Q distribution of the nominal sample in the data. 



TABLE I. Parameters of our fit to the Q distribution 



Parameter 


Description 


To 


Breit Wigner width of Q signal distribution, T(D*^) 


On 


Mean of Q signal distribution 


Ns 


Number of signal events 


fmis 


Fraction of mismeasured signal 




Resolution on measured Q for mismeasured signal 




Number of background events 


k 


cjQ scale factor, fixed to 1 




Coefficients of background polynomial, fixed from simulation 



We have many other parameters of the fit that can be varied or allowed to float for testing 
purposes. We can allow a scale factor on each candidate's aq to model a systematic mistake 
in our tracking system caused, for example, by not properly accounting for the material of 
the detector. In our standard fits we only allow the normalization of the background to float, 
but we can either vary the shape as indicated by the simulation or allow the parameters of 
the background polynomial to float as a measure of the small systematic uncertainty due to 
the background shape. 

Table | summarizes the parameters of our fit. Note that the ctq scale factor k and the 
background shape parameters -81,2,3 are fixed in our nominal fits. We minimize the likelihood 
function 

N 

L = 2{N, + N,)-2J2 log[NsS{Qi, gq,- Tq, Qo, frms. o^,,, k) + N,B{Q,; 3,^2,3)], (15) 

i=l 

where S and B are respectively the signal and background shapes discussed above. 

The fitter has been extensively tested both numerically and with input from our full 
simulation. We find that the fitter performs reliably giving normal distributions for the 
floating parameters and their uncertainties. It also reproduces the input r{D*~^) from to 
130 keV. Its behavior on each of the three data samples: nominal; tracking selected; and 
kinematic selected in the full simulation is discussed below. We note that if all the parameters 
are allowed to vary simultaneously there is strong correlation among the intrinsic width, Fq, 
the fraction of mismeasured events, fmis, and the ctq scale factor, k, as one would expect. 
Thus our nominal fit holds k fixed, but in our systematic studies we either fix one of the 
three or provide a constraint with a contribution to the likelihood if the parameter varies 
from its nominal value. 

IV. FIT RESULTS 

As a preliminary test to fitting the data we run the complete analysis on a fully simulated 
sample that has about ten times the data statistics and is generated with a range of under- 
lying r{D*^) from to 130 keV. We do this for nominal, tracking, and kinematic selected 
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TABLE II. Results of the fits described in the text. The fit parameters are summarized in 
Table |. The uncertainties are statistical. 





Sample 


Parameter 


Nominal 


Tracking 


Kinematic 


To (keV) 


98.9 ± 4.0 


106.0 ± 19.6 


108.1 ±5.9 


Qo (keV) 


5853 ± 2 


5854 ± 10 


5850 ± 4 


Ns 


11207 ± 109 


353 ± 20 


3151 ±57 


fmis (%) 


5.3 ±0.5 


NA 


NA 


amis (keV) 


508 ± 39 


NA 


NA 


Nb 


289 ± 31 


15 ±7 


133 ± 16 



samples. For the nominal sample we note that the fit is not stable if all the parameters are 
left to vary freely. We have found that if we constrain the fraction of mismeasured signal to 
(5.3 ± 0.5)% as indicated by the simulation over the range of generated widths of the D*~^ 
then we get a stable result. This constraint makes the fit to the simulated nominal sample 
have no significant offset between the generated and measured values for the width of the 
D*^. The tracking and kinematic selected samples have a negligible amount of mismeasured 
signal according to the simulation and in fits to these samples we fix fmis to zero. These 
simulated samples are also consistent with no offset between the generated and measured 
values for the width of the D*~^. We also note that in all three simulated samples there 
are no trends in the difference between measured and generated width as a function of the 
generated width; the offset is consistent with zero as a function of the generated width of the 



D*^. Table |ITl| summarizes this simulation study. We will apply these offsets to the fit value 
that we obtain from the data. For the energy release all samples show small shifts, — 7 ± 3 
keV for the nominal, —12 ± 10 keV for the tracking, and —12 ± 5 keV for the kinematic. 

Figures ^, ^, and [1^ respectively display the fit to the nominal, tracking, and kinematic 
selected data samples. The results of the fits are summarized in Table |I[ Correlations 
among the fioating parameters of the fit are negligible. Figure ^ displays the likelihood as 
a function of the width of the D*~^ for the fits to the three data samples. 

The agreement is excellent among the three fits, and when the offsets from Table |ITI| are 
applied we obtain 

Nominal Sample T{D*+) = 96.2 ± 4.0 keV, (16) 
Tracking Selectedr(L)*+) = 104 ± 20 keV, and (17) 
Kinematic Selected T{D*+) = 103.8 ± 5.9 keV. (18) 



The data sample and results are summarized in Table |T|. The uncertainties are only statis- 



tical. We discuss systematic uncertainties in the next section. 
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FIG. 8. Fit to nominal data sample. The different contributions to the fit are shown by 
different colors. 



TABLE III. Summary of the data sample, simulation biases, and fit results. 





Sample 


Parameter 


Nominal 


Tracking 


Kinematic 


Candidates 
Background Fraction (%) 

Tfit — Tgenerated (keV) 

Fit Fo (keV) 
D*+ Width (keV) 


11496 
2.51 ±0.27 
2.7 ±2.1 
98.9 ± 4.0 
96.2 ±4.0 


368 
4.1 ± 1.9 
1.7 ±6.4 
106.0 ± 19.6 
104 ± 20 


3284 
4.05 ± 0.49 
4.3 ±3.1 
108.1 ±5.9 
103.8 ±5.9 
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FIG. 9. Fit to tracking selected data sample. The different contributions to the fit are shown 
by different colors. 
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FIG. 10. Fit to kinematic selected data sample. The different contributions to the fit are 
shown by different colors. 
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FIG. 11. Likelihood function versus measured D*^ width for the nominal (a), tracking (b), 
and kinematic (c) selected data samples. 



V. SYSTEMATIC UNCERTAINTIES 

We discuss the sources of systematic uncertainties on our measurements of the width of 
the D*~^ in the order of their size. The most important contribution is the variation of the 
results as a function of the kinematic parameters of the D*~^ decay as shown most clearly in 
Figure |^. We estimate this uncertainty by repeating the fits described above in three bins 
for each of the six kinematic parameters and taking the uncertainty as the largest observed 
variation from the nominal values in Table 0. We obtain uncertainties of ±16 and ±8 keV 
on T{D*~^) and Qo respectively. 

The next most important contribution comes from any mismodeling of ag's dependence 
on the kinematic parameters. We estimate this by varying our cut on ctq from 75 to 400 from 
the nominal 200 keV and repeating our analysis with all parameters fixed except allowing the 
error scale factor k to vary freely. This indicates that the resolution is correct to ±4%, and 
we then repeat our standard analysis with k fixed at 0.96 and 1.04. We find uncertainties of 
±11, ±9, and ±7 keV on r{D*~^) for the nominal, tracking, and kinematic selected sample. 
For Qo this uncertainty is negligible except in the tracking selected sample where it is ±4 
keV. 

We take into account correlations among the less well measured parameters of the fit, 
such as k, fmis, and amis, by fixing each parameter at ±la from their central fit values, 
repeating the fit, and adding in quadrature the variation in the width of the D*~^ and Qo 
from their central values. We find uncertainties of ±8, ±9, and ±9 keV on the width of the 
D*^ for the nominal, tracking, and kinematic selected sample, and respectively ±3, ±4, and 
±5 keV on Qq. 

We have studied in the simulation the sources of mismeasurement that give rise to the 
resolution on the width of the D*~^ by replacing the measured values with the generated 
values for various kinematic parameters of the decay products. We have then compared these 
uncertainties with analytic expressions for the uncertainties. The only source of resolution 
that we cannot account for in this way is a small distortion of the kinematics of the event 
caused by the algorithm used to reconstruct the origin point described above. This 
contributes an uncertainty ±4 keV on the width of the D*~^ and ±2 keV on Qo- We have 
also checked that our simulation accurately models the line shape of other narrow resonances 
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visible in our data. Notably the decay A° —>■ p7T~, has a Q only seven times that of D*'^ 
-^°^siow ^he A° decay we select the vr" to have a momentum in the range of those in the 
D*~^ decay, and the visible widths agree to a few percent between data and simulation. 

We consider uncertainties from the background shape by allowing the coefficients of the 
background polynomial to float. We observe changes on the width of ±4 keV for the nominal 
sample and ±2 keV for the tracking and kinematic selected samples. We have also released 
our kinematic selection cuts which causes the background to increase by a large factor. This 
causes a change which is small compared to allowing the coefficients of the background shape 
polynomial to float. Variations in the background have a negligible effect on Qq. 

Minor sources of uncertainty are from the width offsets derived from our simulation and 
given in Table |T|, and our digitized data storage format which saves track parameters with 
a resolution of 1 keV and contributes an uncertainty of ±1 keV on the width of the D*~^ and 

Qo. 

An extra and dominant source of uncertainty on Qo is the energy scale of our measure- 
ments. We evaluate this uncertainty by selecting Kg tt+tt^ decays in our data. The 
daughters tracks of the Kg candidates are required to pass the same selection criterion as 
those described above in the nominal sample, the decay vertex is required to be inside the 
beam pipe, and the vertex is required to be signiflcantly separated from the overall event 
vertex. Our Kg sample is quite clean, less than 1% background under the mass peak, and has 
millions of candidates. We then plot the mean of the tt+tt" invariant mass as function of the 
momentum of the daughters. We flnd that above a daughter momentum of 500 MeV/c the 



Kg mass agrees with its expected value |jTT|. We apply corrections, less than 0.3% relative, to 
tracks between 100 and 500 MeV/c to bring the mass peaks into agreement with the nominal 
value. These corrections only affect the slow pion and produce a shift in Qo of —4 keV and 
a negligible change in the width. We evaluate uncertainties in the energy scale by varying 
an overall momentum scale to give a ±30 keV variation, the uncertainty, of the Kg tt+tt^ 
mass, and applying the statistical errors we obtain from the calculations of the momentum 
corrections discussed above. Conservatively we add in quadrature twice the observed shift. 
We observe an uncertainty of 8 keV on Qo and 1 keV on the width due to uncertainty in the 
energy scale of our measurements. 

Table fV\ summarizes the systematic uncertainties on the width of the D*~^ and Qo- 



VI. CONCLUSION 

We have measured the width of the D*~^ by studying the distribution of the energy 
release in D*~^ —>■ D^tt^ followed by — >■ K~'i:^ decay. We have done this in three 
separate samples, one that is minimally selected, a second that reduces poorly measured 
tracks due to misassociated hits and non-Gaussian scatters in the detector material, and 
a third that takes advantage of the kinematics of the decay chain to reduce dependence 
on mismeasurements of kinematic parameters. The resolution on the energy release is well 
modeled by our simulation, with agreement between the sources of the resolution as predicted 
by the simulation and analytic calculations. The largest sources of uncertainty are imperfect 
modeling of the dependence of the mean energy release on the kinematics of the decay chain, 
the simulation of the error on the energy release, and correlations among the parameters of 
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TABLE IV. Systematic uncertainties on the width of the D*~^ and Qq 





Uncertainties in keV 
Sample 




Nominal 


Tracking 


Kinematic 


Source 


6T{D*+) 


6Qo 


6T{D*+) 


6Qo 


6T{D*+) 




Dependence on Kinematics 


16 


8 


16 


8 


16 


8 


Mismodehng of cjq 


11 


< 1 


9 


4 


7 


< 1 


X lb V_y^Ji i ClCl/ljlUllO 


8 


3 


9 


4 


9 


5 


Vertex Reconstruction 


4 


2 


4 


2 


4 


2 


Background Shape 


4 


< 1 


2 


< 1 


2 


< 1 


Offset Correction 


2 


3 


6 


10 


3 


5 


Energy Scale 


1 


8 


1 


8 


1 


8 


Data Digitization 


1 


1 


1 


1 


1 


1 


Quadratic Sum 


22 


12 


22 


16 


20 


14 



the fit to the energy release distribution. With our estimate of the systematic uncertainties 
for each of the three samples being essentially the same we chose to report the result for the 
sample with the smallest statistical uncertainty, the minimally selected sample, and obtain 

r(D*+) = 96±4±22 keV, (19) 

where the first uncertainty is statistical and the second is systematic. We note that if we form 
an average value taking into account the statistical correlations among our three measures 
we get a result that is nearly identical with Equation |l^ since the average is dominated by 
the result with the smallest statistical uncertainty. 

This is the first measurement of the width of the D*~^, and our measurement corresponds 
to a strong coupling 

c/ = 0.59 ±0.01 ±0.07, (20) 

and 

gD'^Dn = 17.9 ±0.3 ±1.9. (21) 

This is consistent with theoretical predictions based on HQET and relativistic quark models, 
but higher than predictions based on QCD sum rules. 

We also measure the mean value for the energy release in D*^ D^n^ decay 

Qo = 5842 ± 2 ± 12 keV, (22) 

where the first error is statistical and second is systematic. Combining this with the mass 
of the charged pion, 139.570 MeV, with an uncertainty less than 1 keV [|lT|, we calculate 

mD*(2oio)+ - = 145.412 ± 0.002 ± 0.012 MeV. (23) 
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This agrees with the value from the Particle Data Group, 145.436 ± 0.016 MeV, from a 
global fit of all flavors of D*-D mass differences. It also agrees well with the best previous 
measure from a single experiment that includes an evaluation of systematic uncertainties 
from ACCMOR at 145.39 ± 0.07 MeV @. 
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